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We propose a superconducting circuit comprising a dc-SQUID with mechanically compliant arm 
embedded in a coplanar microwave cavity that realizes an optomechanical system with a degen¬ 
erate or non-degenerate parametric interaction generated via the dynamical Casimir effect. For 
experimentally feasible parameters, this setup is capable of reaching the single-photon, ultra-strong 
coupling regime, while simultaneously possessing a parametric coupling strength approaching the 
renormalized cavity frequency. This opens up the possibility of observing the interplay between 
these two fundamental nonlinearities at the single-photon level. 


Introduction- Recently, the field of optomechanics 
has undergone a period of rapid advancement, having 
achieved ground state cooling of a mechanical resonator 
[1,2], coherent state transfer between cavity and mechan¬ 
ical oscillators [3, 4], optomechanically induced trans¬ 
parency [5, 6], the generation of squeezed light [7, 8], po¬ 
sition measurement precision close to the standard quan¬ 
tum limit [9-11], and photon back-action evasion [12]. 
However, in all of these experiments a large optomechan¬ 
ical coupling is achieved by the introduction of a strong 
classical drive that effectively linearizes the optomechan¬ 
ical coupling [13]. In contrast, taking advantage of the 
intrinsic optomechanical nonlinearity requires a single¬ 
photon optomechanical coupling strength qq that is larger 
than, or on the order of, the optical cavity decay rate n. If 
in addition is an appreciable fraction of the mechanical 
resonance frequency such that the combined nonlin¬ 
earity parameter g^jnujm ^ 1, then the system is in the 
single-photon ultra-strong coupling regime characterized 
by cavity photon blockade [14, 15]. This regime holds 
promise for achieving the long-standing goal of observ¬ 
ing a macroscopic mechanical oscillator in a non-classical 
state of motion; mechanical steady-states in this regime 
can exhibit sub-Poissonian phonon statistics and negative 
Wigner functions [16-18]. These effects can persist even 
in the presence of thermal fluctuations at experimentally 
accessible temperatures [19]. Recent circuit-QED real¬ 
izations of the optomechanical interaction based on ca¬ 
pacitive [20-22] or magnetic coupling [23] suggest that it 
may be possible to reach this ultra-strong coupling limit 
in the near future. 

In a parallel set of investigations, the interplay between 
the radiation pressure nonlinearity and that of a degen¬ 
erate parametric amplifier (DPA) has been of interest, 
where the inherent squeezing of the cavity mode gives rise 
to a change in cavity photon spectrum to which the op¬ 
tomechanical interaction is sensitive. This squeezing has 
be shown to aid in reaching the strong-coupling regime 
by enhancing the effect of a single-photon in the optical 
cavity [24], and can assist in the cooling [25] and displace¬ 
ment sensitivity [26] of the mechanical resonator, as well 


as in the generation of tripartite entangled states [27]. 
As in other studies, these results focus only on the weak 
coupling regime, and explore the resulting linearized op¬ 
tomechanical Hamiltonian. To date, this approximation 
has not been a limitation, as the only optomechanical 
devices with both intrinsic nonlinear susceptibility and 
radiation pressure coupling are whispering gallery res¬ 
onators (WGR) [27-29], where the optomechanical inter¬ 
action is weak {go/n ^ 10“^); these systems operate well 
outside the single-photon coupling limit. Note that it 
is also possible to generate an effective DPA term from 
the linearized optomechanical Hamiltonian in terms of 
the normal (polariton) modes of the system when using 
a drive detuned to the red-sideband [30-32]. However, 
this transformation necessarily removes the optomechan¬ 
ical interaction. Finding a system that is capable of 
reaching the single-photon, ultra-strong coupling regime 
while simultaneously possessing a degenerate, or non¬ 
degenerate, parametric interaction would allow for inves¬ 
tigating the combined effect of these nonlinearities at the 
single-photon level. 

Motivated by the desire to realize experimentally fea¬ 
sible systems permitting both single-photon ultra-strong 
optomechanical coupling and DPA interactions, in this 
Letter we describe an optomechanical scheme involving 
a dc-SQUID with mechanically compliant arm embedded 
in a coplanar microwave cavity that can in principle real¬ 
ize the single-photon ultra-strong coupling limit, while 
the cavity mode is simultaneously coupled to a DPA. 
This scheme was originally examined in the context of 
optomechanical displacement detection [33, 34], and in 
related subsequent experiments (but without the cavity 
mode) [35, 36]. Here, we take advantage of the recent 
demonstration of the dynamical Casimir effect (DCE) in 
related circuit devices using flux-modulated SQUIDs [37- 
39], and the well-known connection between the DCE in 
a single-mode high-Q cavity and a DPA [40]. In modulat¬ 
ing an external flux applied to the SQUID loop around 
a fixed dc-bias, the sinusoidal variation of the SQUID 
effective inductance in turn modulates the frequency of 
the cavity mode non-adiabatically, leading to photon pro- 
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duction from the quantum vacuum that can be viewed 
as a nonlinear susceptibility. An effective nonlin¬ 
ear medium is generated by amplifying quantum vac¬ 
uum fluctuations in the cavity mode through modulat¬ 
ing its electrical length. Provided that the modulation 
frequency is much higher than the mechanical frequency, 
the resulting separation of timescales allows the optome¬ 
chanical and parametric interactions to decouple, and an 
analysis of experimentally feasible parameters suggests 
that this system is in the single-photon ultra-strong op¬ 
tomechanical coupling regime ^ 1, with an in¬ 

trinsic DPA coupling strength that is on the order of the 
renormalized cavity frequency GHz). 

Hamiltonian- The device, depicted in Fig. 1, consists 
of a coplanar microwave cavity of length I and frequency 
cjc, characterized by its inductance and capacitance per 
unit length Lc and Cc, respectively, bisected by a dc- 
SQUID. The SQUID comprises two Josephson junctions 
with critical current U and capacitance Cj. One arm 
of the SQUID loop is mechanically compliant, forming a 
doubly clamped resonator of length /qsc- We consider 
only the fundamental mechanical mode displacements 
in the plane of the SQUID loop and assume that the 
mechanical resonator can be modeled as a harmonic os¬ 
cillator with the y coordinate giving the center-of-mass 
displacement. The magnetic flux threading the SQUID 
loop is given by ^ext(^) = ^ext(O) + AHextGsc^, where 
^ext(O) = ^ext is the flux with the mechanical oscillator 
fixed at ^ = 0, Hext is the local magnetic field in the vicin¬ 
ity of the mechanical resonator, and A is a dimensionless 
geometrical factor that accounts for the nonuniform dis¬ 
placement of the resonator along its extension. 

To model the SQUID, we take the sum and difference of 
the gauge-invariant phases, and 02, across the Joseph¬ 
son junctions, 7 ± = (0i ± 02 )/ 2 , and use the phase-field 



FIG. 1. (color online) Layout of a dc-SQUID with a me¬ 
chanically compliant segment (red) embedded in a coplanar 
microwave cavity of length /, and driven by a time-dependent 
external flux bias Text(t) (orange). The SQUID, here greatly 
enlarged with respect to the cavity for visibility, is assumed to 
be a lumped element at x = 0, with two identical Josephson 
junctions (blue) characterized by their critical current U and 
capacitance Cj. 


coordinate 0(x, t) for the microwave cavity. The current 
and voltage along the cavity are given by the telegraph 
relations. 


V{x,t) 


$0 d(f>{x,t) 
27rLc dx 
^0 50(x,t) 

27r dt 


(1) 

( 2 ) 


where <Fo = hj(2e) is the flux quantum. Assuming that 
the SQUID can be modeled as a lumped circuit element 
at the center x = 0 of the cavity, the closed system equa¬ 
tions of motion for the cavity, SQUID, and mechanical 
oscillator are given by [33] 

dt^ LcCc dx‘2 ’ ^ ’ 


cos ( 7 +)sin( 7 _) 


+ 2^, 


-1 


7_ — TT n + 


^ext T AHgxt^ oscV 


^ 7 + + sin ( 7 +) cos ( 7 -) + 


$0 

$0 d(j)(Q,t) 

AttLcIc dx 


= 0 , 


and 


( 4 ) 
= 0 , 

( 5 ) 


my + - ;^ABext^osc7- = 0 , ( 6 ) 


where wj = 727r/c/ (Cj^o) is the Josephson junction 
plasma frequency, Pl = 27rL/c/^o is a dimensionless 
parameter characterizing the SQUID self-inductance L, 
and the integer n arises from the requirement that the 
phase around the SQUID loop must be single valued. 
Equation (3) is the wave equation for the microwave cav¬ 
ity, Eq. (4) describes the current circulating through the 
SQUID loop, Eq. (5) gives the average current in the 
loop, and Eq. ( 6 ) is Newton’s equation for the mechan¬ 
ical resonator driven by the Lorentz force. The current 
and voltage across the SQUID must also obey the bound¬ 
ary conditions 

d(p{±l/2,t) ^ d(p{0-,t) ^ d(p{0+,t) 

dx ’ dx dx ^ ^ 

and 

. L d'^(f>{0,t) _ d4>{0~,t) d(f>{0+,t) 

41/c dtdx dt dt ’ 

respectively. 

Using Eqs. (3)-( 8 ), we can derive approximate equa¬ 
tions of motion that govern the interaction between the 
cavity and mechanical resonators that are determined 
solely by the SQUID parameters and boundary condi¬ 
tions. In what follows, we assume that the SQUID 
plasma frequency ujj satisfies ujj ujc ^ cjm, allow¬ 
ing us to ignore the SQUID inertia terms in Eqs. (4) and 
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(5) . Assuming that ^ 1 (i.e. the self-inductance of 
the SQUID loop is negligible), together with dropping the 
SQUID inertia terms, allows us to eliminate the SQUID 
phase coordinates 7 ± from the equations of motion. Fur¬ 
thermore, we assume that the cavity current is much less 
than the Josephson critical current, \I/Ic\ ^ 1, and that 
the change in flux due to the small displacements of the 
mechanical resonator is much less than the flux quantum 
|A5ext^osc^/^o| ^ 1- Keeping terms to first order in 
and to leading second order in /, the equation of motion 

(6) for the mechanical oscillator becomes approximately 


my -h muJmV = 


A^ext^o 


4/c 


■ sec 


7r<F 


ext 




tan 


TT^e 




( 9 ) 

and the voltage boundary condition (8) can be expressed 
as 




^0 ' d<f){Q ,t) _ d(j){0+,t) 
27r dt dt 


( 10 ) 


with 


U (4>ext5 y) — Lj (4>ext) 1 


A-Bext^osc^ , / TT^ext \ 


(11) 

where we have kept only the current independent contri¬ 
bution to the SQUID effective inductance [the next term 
being of order (///c)^] 


5^ cos{uJdt). To first order in the small oscillation am¬ 
plitude the effective length parameter becomes 


T (4>ext5 


Lj (^dc) 

Lc 


1 + 


—— tan 
4>o 


f 7r4>dc 

V ^0 


'^-^ext^osc^ 

(^oA) 


tan 



cos (cddQ 


(16) 


Restricting ourselves to a single cavity mode (the exten¬ 
sion to multiple cavity modes is straightforward), and 
temporarily fixing the position of the mechanical oscilla¬ 
tor at ^ = 0, then the cavity phase field that satisfies the 
current boundary conditions (7) is 


4i{x,t) 


COS [ko{t) {x1/2)] X < 0, 
+(/)(t) cos [ko{t) {x — 1/2)] X > 0, 


(17) 


with the wavenumber ko{t) determined by the flux 
boundary condition (15): 



LJ 

Tj(^dc) 



<I>o 


tan 


f 7r<Fdc 

V ^0 


cos (uJdt) 


(18) 

For LJ/Lj{^dc) ^1/2, the left hand side of Eq. (18) is 
well approximated by {kJ/2)‘^ that, together with uj‘1 = 
k^/[LcCc) from Eq. (3), allows the cavity mode frequency 
to be expressed as 


L j{^ ext) 


4>q 

47r/c 


sec 


/ TT^ext \ 


( 12 ) 


^l{t) = (71 


7r^<F 

—— tan 
<Fo 


f 7r4>dc 

V ^0 


cos (cddQ 


, ( 19 ) 


The validity of the approximations used in the derivation 
of the system equations of motion requires not approach¬ 
ing too close to the half-integer flux bias point: 



/3l sec 


f TT^ext \ 

V ^0 y 

/ TT^ext \ 

v^/ 


«1, 

< 1 . 


(13) 

(14) 


Taking the time integral of both sides of (10), setting 
the integration constant to zero, and inserting the cavity 
current expression (1), Eq. (10) for the phase across the 
SQUID can be written as a Robin-type boundary condi¬ 
tion: 


l^(^ext,l/) d(p{0,t) 
Lc dx 


(f>{o+,t) - 4>{o ,t)- 


(15) 


where (cdc*^)^ = 4/[Cc/Tj(^dc)]- This expression is sim¬ 
ilar to that of Ref. [39], where this flux-dependent fre¬ 
quency modulation was used in generating radiation via 
the DCE in a dc-SQUID array. The dependence of 
on the SQUID inductance (11) evaluated at the 
dc-flux bias gives a cavity frequency renormalization 
(X cos (7r^dc/4>o)- The plasma frequency ujj is also renor¬ 
malized by an amount proportional to a/cos (7r$dc/4>o), 
and therefore the plasma frequency remains well above 
for all flux biases compatible with Eqs. (13) and (14). 

Substitution of the phase field (17) into the oscillator 
equation of motion (6) gives the cavity force acting on the 
resonator when ^ = 0. Under our assumptions of small 
and slow mechanical displacements, the Lorentz force in 
(9) is unchanged to good approximation and the equation 
of motion for the mechanical oscillator becomes 


In the absence of a mechanical oscillator, this bound¬ 
ary condition is of the same form as that used in super¬ 
conducting realizations of the DCE, where the quantity 
L (4>ext7 0) /Lc can be viewed as a flux-tunable length pa¬ 
rameter [37, 41]. 

Let us now suppose that the external flux is weakly 
modulated around some fixed dc bias $ext = 4>dc + 


my + mLol,y J[^) CJsiv? (20) 

^ A^ext^osc Lj (i»dc) 

(i>o/27r) Ld 

X tan (71^11, 










































4 


where we have assumed ujd ^ ujrn such that terms pro¬ 
portional to cos{uJdt) average to zero over a single me¬ 
chanical oscillation; modulation of the applied flux does 
not significantly affect the coupling between the cavity 
and mechanical resonator. Equation (20) allows us to 
determine expressions for both the mechanical portion 
of the system Lagrangian and that of the interaction 
with the cavity. The remaining cavity terms in the La¬ 
grangian follow from the wave equation (3), and the total 
Lagrangian for the system can be written as 

£ (4>, y, 4>, y) = (21) 

1 '^-^ext^osc^ Lj (^dc) 

^2 i^o/^TT) LJ 

X tan 

where the effective phase mass is defined to be 



Defining the lowering and raising operators for the cav¬ 
ity mode (d,d+) with respect to as well as those 
for the mechanical resonator (6, 6+) satisfying the usual 
bosonic commutation relations, the system Hamiltonian 
takes the form 




H = 

ha 


a -h hujrnb'^b 


( 22 ) 


- — -h (a -h a+)^ 

-^{a + a+)jb + bj 
+ HE (ae+*‘^>’* + 0+6“*“^’*) , 


where, for completeness, we have included a term corre¬ 
sponding to the pumping of the cavity by a classical mi¬ 
crowave field at frequency ujp with amplitude E. Here, 
the coupling strengths are expressed as 


a = 


cjf 7r5^ 
4 <Lo 


tan 



(23) 


and 


go = , tan —— , 

" (^oA) LJ V ^0 y 


(24) 


where Uzp = yh/JjmuJrn) is the zero-point displace¬ 
ment of the mechanical resonator. The coupling a is the 
product of both the flux-modulation amplitude [e.g. see 
Eq. (19)] and the intrinsic parametric coupling y. Eac- 
toring out the former, allows us to define the effective 
parametric coupling strength resulting from the DCE to 
be X = 


Moving to a frame rotating at uJd/^^ and dropping off- 
resonant terms we obtain 


H = —hAa'^a -f hujmb'^b — 


hgoa'^a (h + 6^^ + HE (de 


a 

—iSt 


+ 

+ a+e+®'’*) 


(25) 


with detunings A = and 6 = Ud/^—Up. Equa¬ 

tion (25) describes the radiation-pressure interaction be- 
tween mechanical and microwave oscillators driven by an 
effective DPA term arising from the DCE, as well as by a 
classical linear pumping term. In addition to providing a 
resource for cavity mode squeezing, the DPA can also be 
utilized to control the mechanical resonator. By detun¬ 
ing the driving frequency by an amount proportional to 
half the mechanical resonance frequency, A = ±ci;^/ 2 , it 
becomes energetically favorable for a pair of photons gen¬ 
erated by the DCE to enter the cavity by emitting (ab¬ 
sorbing) a phonon into (from) the resonator, thus heat¬ 
ing (cooling) the mechanical mode. Similarly, detuning 
by A = EuJrn makes this emission or absorption a two 
phonon process. 

Although we have restricted ourselves to a single cavity 
mode, the extension to multiple modes is straightforward. 
In the case of two cavity modes, modulating the cavity 
resonance frequency at the sum of the two dc-biased fre¬ 
quencies gives rise to a resonant non-degenerate para¬ 
metric coupling (dc + d+c+), where the operators (c, c+) 
correspond to a cavity mode with non-vanishing current 
at the location of the SQUID, i.e. the cavity mode phase 
must satisfy 7 ^ 0 at x = 0. In addition, driv¬ 

ing at the difference of the mode frequencies results in 
a beam-splitter interaction (dc+ + d+c), which has also 
been investigated for membrane in the middle configura¬ 
tions [42, 43]. With the mechanical resonator coupled to 
both modes via radiation pressure coupling, possibly in 
the single-photon strong-coupling limit, such a configu¬ 
ration provides a rich testbed for investigating tripartite 
entanglement in strongly nonlinear systems [44, 45]. 

Ultra-strong coupling regime- Given that the optome¬ 
chanical coupling is dependent on the external magnetic 
field, it is advantageous to consider a device fabricated 
using niobium, with critical field Be ^ 198 mT, rather 
than aluminum {Be ^ 10 mT). A set of experimen¬ 
tally feasible parameters compatible with this choice is: 
4>dc = 0.354>o, cjy^(O) = 27rX 10 GHz, = 27rx 10 MHz, 
m = 10“^^ kg, lose = 10 /im, A = 2 / 7 r, Hgxt = 40 mT, 
le = 100 nA and Lei = I nH. At this dc-bias, the 
SQUID renormalizes the cavity frequency to nearly half 
of its original value 2 tt x 4.5 GHz. In addition, 

we assume quality factors for the dc-biased cavity and 
mechanical resonator to be = 5 x 10 ^ and Qm = 10 ^, 
respectively. The former value has already been achieved 
in flux-tunable A1 resonators [46], although the addition 
of a flux bias degrades the quality factor in these oscil¬ 
lators. In contrast, the quality factor for Nb cavities is 
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System 

n 


90 

LOm 

9o 

KiLU'jji 

A 

27T 

Microwave LC-drum [4] 

6- 

lO-^ 

2-10“® 

10“® 

n/a 

Si zipper cavity [7] 

2 • 


3-10“^ 

6•10"® 

n/a 

Circuit-QED qubit [22] 

4- 

10“^ 

2 • 10“^ 

8 • 10“'‘ 

n/a 

cCPT - resonator [20] 

10 


1 

10 

n/a 

Stripline - cantilever [23] 

20 

a 

- 

- 

n/a 

Si3N4 WGR [27] 

2 • 

10“^ 

2•10“^ 

4•10“^ 

110 Hz 

SQUID - resonator 

13 


0.1 

1.3 

1 GHz 


^ Ref. [23] assumes a cavity quality factor of 10® taken from 2>D 
superconducting cavity geometries rather than the stripline 
configuration considered in that work. As such, the realizable 
values for qq are likely to be lower. 

TABLE I. Ratio of the optomechanical coupling strength 
to that of the cavity decay rate, go/n and mechanical fre¬ 
quency combined quantum nonlinearity parameter 

g^jKUJra^ and intrinsic parametric coupling strength x for a 
selection of recent experimental (above dashed-line) and the¬ 
oretical (below dashed-line) optomechanical systems incorpo¬ 
rating macroscopic mechanical resonators. 


marginally affected by an external bias provided that the 
renormalized frequency is large enough for thermal effects 
to be ignored [47]. For the frequencies considered here, 
the thermal occupation of the cavity mode is ^ 10“^^ 
at 10 mK, and thus we expect a negligible change in 
the quality factor. Substitution of these parameters into 
Eq. (24) gives g^j^ ~ 13 and g^jKUJjn ~ 1.3, indicating 
that this setup lies within the single-photon ultra-strong 
coupling regime. In addition, the intrinsic parametric 
coupling strength at this dc-bias is y = 27r x 1 GHz. The 
total parametric coupling strength a (23) will be much 
lower due to the weak modulation amplitude (54>/4>o ^ 1 
and the stability requirements [25, 48] of the system. A 
comparison of these figures of merit with recent experi¬ 
mental and theoretical optomechanical systems incorpo¬ 
rating macroscopic mechanical resonators is presented in 
Table I. 

Although we have tacitly assumed that our optome¬ 
chanical system involves only a single mode of the me¬ 
chanical oscillator, in practice all odd harmonics that are 
symmetric about the midpoint of the resonator will cou¬ 
ple to the SQUID through the modulation of the exter¬ 
nal flux. Our single-mode approximation rests on the as¬ 
sumption that higher-order modes are spectrally well sep¬ 
arated from the fundamental. However, given the large 
single-photon coupling strength in this device, it is of in¬ 
terest to ask whether additional higher-modes are also 
strongly coupled to the cavity. The tripartite interac¬ 
tion between a cavity and two mechanical modes has al¬ 
ready been explored for generating entanglement [49, 50], 
two-mode squeezing [51], and EPR states [52], while as 
an ensemble of mechanical modes has been proposed for 


quantum information processing [53] . Assuming that the 
higher modes are harmonics of the fundamental, the ra¬ 
diation pressure coupling (24) is reduced by a factor of 
l/^3/2, ^ labels the harmonic, arising from a re¬ 

duction in both the geometrical factor accounting for the 
net change in SQUID loop area and a decrease in the 
zero-point amplitude at higher frequencies (assuming the 
mass of each mode is approximately the same). For the 
parameters detailed here, this estimation suggests that 
up to the first six odd-harmonics of the resonator can be 
simultaneously strongly (but not ultra-strongly) coupled 
to the microwave cavity. 

Conclusion- We have shown that a dc-SQUID with 
mechanically compliant segment embedded in a copla- 
nar microwave cavity modulated by a time-dependent 
external flux gives rise to an optomechanical interaction 
where the cavity mode is coupled to a degenerate or non¬ 
degenerate parametric amplifier. Using experimentally 
feasible device parameters indicates that this setup is ca¬ 
pable of operating in the single-photon ultra-strong op¬ 
tomechanical coupling regime, while simultaneously pos¬ 
sessing a degenerate parametric coupling strength on the 
order of the dc-bias renormalized cavity frequency. Addi¬ 
tional higher-modes of the mechanical resonator can be 
strongly-coupled simultaneously in this setup. In con¬ 
trast to previous proposals for strongly-coupled optome- 
chanics, where the frequency of the mechanical resonator 
is intrinsically low [23], or is renormalized downward 
[20, 21] such that thermal effects conspire to degrade the 
quantum signatures in these setups, the present scheme 
lies within the parameter space where quantum signa¬ 
tures such as negative Wigner functions are predicted 
to persist at standard dilution refrigerator temperatures 
[19], and thus is well-suited for investigating macroscopic 
quantum states of the mechanical resonator. Compared 
to current WGR systems, the device presented here gives 
a four and nine order of magnitude increase in go/n and 
gQlnujm^ respectively. Combined with the strong intrin¬ 
sic parametric coupling x, this setup opens the door to 
exploring the expected rich physics arising from the in¬ 
teraction between these two principle nonlinearities, in 
a single or multi-mode configuration, at the individual 
photon level. 
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